Variety of c-axis collective excitations in layered multigap superconductors 
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We present a dynamical theory for the phase differences along a stacked direction of intrinsic 
Josephson junctions (UJ's) in layered multigap superconductors, motivated by the discovery of 
highly-anisotropic iron-based superconductors with thick perovskite-type blocking layers. The dy- 
namical equations describing AC and DC intrinsic Josephson effects peculiar to multigap UJ's are 
derived, and collective Leggett mode excitations in addition to the Josephson plasma established 
in single-gap UJ's are predicted. The dispersion relations of their collective modes are explicitly 
displayed, and the remarkable peculiarity of the Leggett mode is demonstrated. 

PACS numbers: 74.50.+r,74.70.Xa 



Highly-anisotropic layered High-T c superconductors 
are natural nano-scale stacks of Josephson junctions, i.e., 
intrinsic Josephson junction (IJ J) arrays, since supercon- 
ducting and insulating layers with atomic thickness reg- 
ularly alternate along the crystalline c-axis. Their high- 
quality single crystals clearly exhibit Josephson effects 
only in c-axis electromagnetic response, which are called 
intrinsic Josephson effects (IJE's). IJE's have been ex- 
perimentally confirmed in various layered High-T c copper 
oxide superconductors, such as Bi 2 Sr 2 CaCu208 (lUsl- 
An intriguing feature in IJE's is unique dynamics aris- 
ing from couplings between the stacked junctions. Two 
types of inter- Junction couplings due to inductive [3, S| 
and capacitive 0, 0] origins have been mainly proposed. 
The inductive coupling in UJ's is very strong since the 
in-plane magnetic penetration depth characterizing the 
magnetic field screening range extends over several hun- 
dred junctions. Meanwhile, the capacitive one is not so 
strong since the charge screening length is comparable to 
the layer thickness. However, it has a significant role on 
IJE's due to the atomic-scale structure [8|. 

Most of High-r c cuprate materials are identified as 
single-band superconductor. One then defines just a 
single phase difference between consecutive supercon- 
ducting layers. The dynamics of the phase differ- 
ence in single-gap UJ's has been intensively studied 
after the discovery of High-T c cuprate UJ's Q. In 
this paper, we extend the dynamical theory for the 
phase difference to multi-gap UJ's, in which more 
than one phase differences are active through stacked 
all junctions. Our motivation comes from the dis- 
covery of highly-anisotropic iron-based supercondcu- 
tors, such as (Fe2As2)(Sr4V 2 06), (F e 2 P2)(Sr 4 Sc206) and 
(Fe 2 As2)(Sr 4 (Mg, Ti) 2 6 ) [TMl. These compounds 
contain thick perovskite-type blocking layers Sr4M2 06 
(M = Sc, Cr, V) with a thickness of ~ 15 A, which 
clearly remind us of Bi 2 Sr 2 CaCu208- The first princi- 



ples calculations on these materials indicate that they are 
multiband systems with strong two-dimensional charac- 
ter [l3], EH whose anisotropy is estimated to be compa- 
rable to Bi 2 Sr 2 CaCu20 8 . In fact, the experimentally ob- 
served resistivity of their polycrystalline samples [10j ex- 
hibits a broad superconducting transition in the presence 
of external magnetic field, which is a clear sign of high 
anisotropy. We also note a direct report that single crys- 
tals of an iron-based superconductor, PrFeAsOo.7, show 
the I-V characteristics peculiar to Josephson junctions 
in the c-axis [lj| . 

What is the most fundamental issue in multi-gap UJ's? 
Since the tunneling channel is also multiple, the number 
of collective modes in the phase oscillation is simply ex- 
pected to be multiplied. Confining ourselves to the sim- 
plest two-gap systems, we study the multiple collective 
modes. First, we derive coupled equations of motion for 
the phase differences describing AC and DC multigap 
IJE's. Then, a mode analysis on them clarifies that the 
in-phase mode corresponding to the Josephson plasma is 
not significantly altered while the out-of-phase one, i.e., 
Leggett mode suggested in the presence of two superfluid 
orders by Leggett [16(, emerges as a unique mode. An in- 
triguing focus in this paper is that the c-axis Leggett 
mode is weakly dispersive and favors synchronous oscil- 
lations along c-axis. Such a behavior is striking contrast 
to the Josephson plasma. 

Consider the two-gap UJ's composed of N junctions 
as shown in Fig. [Ha). We assume the pairing interaction 

as - Z)ij=i,2 dijWu^ui^ijl^m on eacn superconduct- 
ing layer, where ipHcr is the electron field operator with 
spin a in the ith band on the Ith superconducting layer. 
The coupling constants 511 and (722 (312 = 521) denote 
the intra-band (inter-band) pairing interaction constants. 
The inter-band pairing interaction generates the inter- 
band Josephosn coupling energy, v ln cos(ip]p —Jf^), in 
the effective action of superconducting phases [17j|. Here, 
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FIG. 1: (color online) (a) Schematic diagram for IJJ's with 
multiple gap superconducting layers, (b) Double-sided type 
junction and the associated boundary condition for ipi. (c) 
Mesa type junction and the associated boundary condition 
for tpi. 



(pi is the phase of the superconducting gap in the ith 
band on the ith superconducting layer, and the coupling 
constant V[ n is given as V[ n = 4Ki n |<7i2/g||AW||A( 2 )|, 
where |A^| is the zth superconducting gap amplitude 
and g = 311322 — (<?12) 2 - The coefficient K; n is the sign 
factor defined clS /^in — 

1 for gi2 > and K in = -1 
for gi2 < 0. For the Josephson coupling between con- 
secutive superconducting layers, one can derive the so- 
called Josephson coupling energy on two channels as 

Ei=i, 2 (^M/ e *) C0S [^+i - V/ W " {e*d/hc)A z l+l l ], where 
j Cj i and j Cj 2 are the Josephson critical current densities, 

j4f +1 1 = (1/rf) fl"! A z (z)dz is the z-component of the 
vector potential, and e* = 2e. Here, we neglect the inter- 
band crossing channel because it is the forth-order pro- 
cess in the coherent tunneling case. On the basis of this 
result, we propose an effective Lagrangian for the two- 
gap IJJ's as 
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where g M = {h/e*)d t <pf ] + A?, v M = {hc/e*)d xV f ] - Af , 
(e*d/hc)A z +1 j, J in = e*Vi n s/h, 
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and ^ = i/Jj — tp\ z> . The parameters, s and d, arc the 
thicknesses of the superconducting and insulating layers, 
respectively, fa (\ a b,i) is the charge screening length (in- 
plane penetration depth) relevant to the ith-band elec- 
trons, e is the dielectric constant in the insulating lay- 
ers, and A° and Af are, respectively, the scalar poten- 
tial and the x-component of the vector potential on the 
Ith superconducting layer. Without losing generality, we 



consider only the z-component of the electric field and 
only the y-component of the magnetic field, which are 
expressed as Ef +ll = -(l/c)d t Af +1<l - {A° +1 - A°)/d 
and Bf +1 1 = (Af +1 —Af)/d— d x A* +1 v As in the single- 
gap IJJ's one can define the inductive [1, HJ and capac- 
itive [(3, 0] coupling constants in the dimensionless form 
as rji = X 2 ab i/sd and a« = e/if/sd for each channel in this 
system. The effective action (pj describes low energy 
dynamics of the superconducting phases in the two-gap 
IJJ's. In the derivation of Eq. (JTJ , all junction parameters 
(e.g. j Ct i and A a b,i) are approximated as local quantities 
for brevity, although they are originally nonlocal ones. 
As for the nonlocal electromagnetic effects, see Ref. fl8j |. 

Now, let us derive the coupled equations of motion 
of the superconducting phase differences. First, one ob- 
tains the so-called Josephson relations associated with 
time and spatial variations of the superconducting phase 
differences as, 

- = ^(1 - aAW)E? +u , (2a) 

d x 6i +1 ,i - = 5^(1 - V& {2) )B? +1 „(2b) 

where a and fj are, respectively, the reduced capacitive 
and inductive coupling constants given as a -1 = a^ 1 + 
a^ 1 and rj^ 1 = 77^ 1 + Tfc , the second-order finite 
difference defined as A^ 2 ' /; + i,; = f 1+2,1+1-1 f 1+1,1+ f 1,1-1 
(for ^ f 1+1,1), <J>o(= 2irhc/e*) the unit flux, and £ = (ai - 
012) /{oil + a 2 ) and £ = (771 — ^/{Vi ~ Vi)- Here, we in- 
troduce 9 l+u = [0i+n+0i+ij)/2 and ipi+ij = tpi+i-tpi. 
Equations (|2a|) and (|2bp are interpreted as the general- 
ized Josephson relations in the two-gap IJJ's. These re- 
lations are reduced to the conventional ones in the single- 
gap IJJ's when £ = £ = 0. On the variation of Af +1 1 , 
the Lagrangian (JTJ) derives the Maxwell equation as 



d * B ?+i,i - - d t E l+i,l = — Ql+l,: 
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and 3?+i,i = 



Jl+1,1 — J2i=l jc,i sm "i + l,Z 

Si=ii{^L j • Here, we add the quasi-particle tunneling 
current j'S^ l , which can be derived microscopically Q . 
Furthermore, we have the continuity equations, which 
can be derived by the variation with respect to ipf 1 0]. 
From the continuity equations with Eq. ([2]) we also have 
the "pseudo" Maxwell equation, which describe the mo- 
tion of the relative phase differences V'z+u, as 

~ e ~ Air 

d * B i+i,i ~ ~ d t E i+i,i = — 2Jin(sini/i ;+ i - smipt) 

Air 
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-Jl+il The "pseudo" electromagnetic fields 
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Equations ([2]), ©, and ([?]) provide a set of equations 
of motion for the phase differences and the electromag- 
netic fields in the two-gap IJJ's, that is, the DC and 
AC Josephson effects in the two-gap IJJ's can be de- 
scribed by these coupled equations. To solve these 
equations it is convenient to use the relation defined as 
ipi = £ m =i ipm,m-i + V'O) where the value of ipo is spec- 
ified as the boundary condition [Figs.fljb) and[]Tc)]. 

Let us focus on the collective phase oscillation modes 
in the two-gap IJJ's. Consider the N junction system 
under the periodic boundary condition along the c-axis. 
For simplicity, the case of £ = C = and j Cj i — j Cj 2 is 
examined in the following. More general cases will be 
published elsewhere. Assuming small oscillations, we lin- 
earize Eqs. Q and ((4]) around 9i+i,i = and ipi = with 



neglecting the dissipation currents ( and l for 
the standard mode analysis [19] . The dynamical simula- 
tion taking into account the quasiparticle contributions 
was performed in Ref. [20| . Eliminating the electric and 
magnetic fields from the coupled linearized equations, we 
can derive the decoupled equations for di+ij and 



(a) 




(b) 




0.65 
0.6 



/ = 5 



0.55 ^ ■ ^ * 

0.5 » — • 
0.45 P"""* 



*'l = 3 

-•1 = 2 
*'] = 1 



(c) 



I) 



<-,A, 



FIG. 2: (color online) (a) Dispersion relations for the 
Josephson-plasma (solid lines) and the Leggett's modes (dash 
lines) when N — 5 and j- ln = 1.0. We set ct\ = a?2 = 0.1, 
771 = 772 = 10 3 , and j c ,i = jc,2- Enlarged views of the dis- 
persion relations at small in-plane wavenumbers for the five 
eigenmodes of the Josephson-plasma (b) and the Leggett's 
modes (c). 
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where A" 2 = A c 2 + A c j, A c 2 = 4j:e*dj c ^/hc 2 , X-J = 

4ire*d\J in \/Hc 2 , 9 = '(02,1, #3,2, • • • , 0n-i,n), and tp = 
'(^2,1, "03,2, iPn-i,n)- The coeffcients, C and L, are 
N x N matrices given as C — (1 + 2a)I — aS and 
L = (1 + 2*7)1 — f}S, where S is an N x N tridiagonal 
matrix with the elements, Su = and j±i = 1, and 
I is the N x N unit matrix. We note that the matrices 
C and L represent, respectively, the capacitive and in- 
ductive couplings between junctions, which are the same 
as those in the single-gap IJJ's. Thus, the collective mo- 
tion of the mean phase differences, which is described 
by Eq. (f5a|) , is understood to be the Josephson plasma. 
Moreover, it is clearly found that its dispersion is brought 
about by the inductive and capacitive couplings between 
junctions. On the other hand, due to two-gap IJJ's, we 
have another collective oscillation mode in the relative 
phase channel, which is described by Eq. (I5bp . In the 
new mode, its origin, i.e., the coupling between junctions 
is found to be induced by the off-diagonal components of 
the matrix N = (1 + 2v)I - vS with 
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We note that the coupling constant 7i n (or v) depends on 
not the inductive and capacitive coupling constants but 
just the inter-band Josephson coupling Jj n . Thus, this 
mode has its origin only in the inter-band pairing inter- 
action. Hence, one understands that Eq. (|5b[) describes 
the Leggett mode in the two-gap IJJ's. From these re- 
sults, it is concluded that the Josephson plasma mode 
is originated from the inductive and capacitive coupling 
arising from the electromagnetic field screening, while the 
Leggett mode is brought about by the intra-layer inter- 
band coupling. 

The dispersion relations of these two eigen-modes are 
obtained from Eqs. (|5a|) and (|5bj) . which are specified in 
terms of the wave numbers k x (in-plane direction) and 
k z = lir/(N + l)d (c-direction) as 



wp(kx,l) = w P (0,/)Wl + 
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with wp(0, 1) = Wpi-y/i + 2a(l — Sf) for the longitudinal 
Josephson plasma and 
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FIG. 3: (color online) (a) Dispersion relations for the 
Josephson-plasma (solid lines) and the Leggett's modes (dash 
lines) when 7i n = 3.0. The other parameter values are the 
same as in Fig. [2] Enlarged viewe at small wave numbers are 
shown in (b) and (c). 



with Al ( 
where si 



2^/v\ c /y/r)i + 772 and wl(0,Z) 



2i/(1 — s{) for the longitudinal Leggett mode, 
= cos{lir/(N + 1)] and u p i and WL Cg are, 
respectively, the Josephson plasma and the Leggett 
mode frequencies, i.e., w p \ — c/yfe\ c and uihcg = 
C\Ja.\ + 02/ -y/eAin • Here, it is clearly found that the ori- 
gin of the Leggett mode is a fluctuation between two 
superfluids which is essential to neutral multiple super- 
fluids. 

We plot the dispersion relations of these eigen modes 
in the case of N = 5 with j Cj1 = j C) 2 for 7j n = 1.0 and 
3.0, respectively, in Figs. [3] and [31 The values of the in- 
ductive and capacitive coupling constants are chosen as 
011 = «2 =0.1 and rjx = 772 = 10 3 , which are the val- 
ues applicable to the cuprate IJJ's. If the Leggett mode 
is a low-energy excitation mode and can lie inside the 
energy gaps as the Josephon plasma, then it is possi- 
ble that both modes are closely located in the low en- 
ergy range as shown in these figures. It should be also 
noted that the Josephson plasma mode with the largest 
c-axis wave number, i.e., I = 5, is the lowest energy one 
close to k x = 0, but this mode changes to the highest 
one for larger values of k x . This is because the large 
inductive coupling, which is predominant in a wide k x 
range, favors 7r phase shift in the phase differences be- 
tween consecutive junctions [2l|. This discussion clearly 
leads to that 7r anit-phase synchronization is preferable 
in the Josephson plasma mode under the presence of the 
layer parallel magnetic field. In fact, strong synchronous 
electromagnetic-wave emission is observed only at the 
zero and weak field in layered High-T c copper oxide su- 
perconductors [22j . On the other hand, the dispersions 
in the Leggett mode does not show such level crossing 
as seen in the figures, since the excitation mode is as- 



sociated with only the density channel. This indicates 
that the Leggett excitation always prefers synchronous 
oscillations along junction stacked direction even in the 
presence of the magnetic field. If the Leggett mode is 
excited by the charge injection or other ways, then the 
synchronized Leggett oscillation emerges and a conver- 
sion into the synchronized Josephson plasma excitation 
due to inherent nonlinearity may occur. 

Finally, we mention that when the difference between 
the two tunneling channels exist (i.e., j c ,i ^ jc.2, ot\ 7^ 
«2, and 771 ^ 772) a mode coupling between the Josephson- 
plasma and the Leggett modes can occur. Such a cou- 
pling effect is an interesting future task. 

In summary, we derived the coupled dynamical equa- 
tions for the phase differences which can be utilized for 
the analysis of AC and DC Josephson effects in the multi- 
gap IJJ's. The equations revealed that multi-gap IJJ's 
have two collective phase oscillation modes, the Joseph- 
son plasma and the Leggett mode whose origins are dif- 
ferent. Moreover, it is found the Josephson plasma and 
Leggett modes favor 7r anti-phase and in-phase synchro- 
nization along the junction stacking , respectively, in a 
wide wave-number range. 
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